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We consider the impact of quantum effects on electron dynamics in a plane linearly polarized 
standing wave with relativistic amplitudes. For this purpose analysis the Lyapunov characteristic 
exponent spectrum with and without allowance for the classic radiation reaction force has been 
analyzed. Based on this analysis it is concluded that the contraction effect of phase space in 
the stochastic regime due to the radiation reaction force in the classical form doesn’t occur when 
the quantum nature of hard photon emission is taken into account. It is shown that electron 
bunch kinetics has a diffusion solution rather than the d’Alambert type solution as in the classic 
description.lt is also revealed that the electron motion can be described using the Markov chain 
formalism. This method gives exact characteristics of electron bunch evolution, such as motion of 
the center of mass and electron bunch dimensions. 

PACS numbers: 42.60.-m 


INTRODUCTION 


Nowadays, much attention is paid to the interaction 
of relativistically strong laser fields with matter, when 
radiation losses play an important role. The interest 
in this subject is connected not only with the funda¬ 
mental problem of charged particle dynamics in elec¬ 
tromagnetic fields, but also with possible applications, 
such as laser-based electron accelerators and hard photon 
sources 0 , 0 . it was recently found that the dynamics 
of charged particles in ultrarelativistic electromagnetic 
fields changes drastically when radiative reaction force is 
taken into account. As was shown in Refs. 01 , at super 
strong laser fields charged particles are drawn into a high- 
field region, thus leading to efficient generation of gamma 
rays and particle beams. The need to consider the back 
impact of the emitted photon !fields! is dictated by the 
facts that the recoil force becomes comparable with the 
Lorentz force and a particle can lose a substantial part 
of its energy in one act of emission. Thus, the effect of 
radiation losses can dramatically change the trajectory 
of the particle over the laser field period. In some cases, 
small changes of particle trajectories over the laser pe¬ 
riod caused by weak radiative losses can be accumulated 
and significantly change particle motion for a time much 
longer than the optical cycle. This follows from the theo¬ 
rem of phase space contraction in the presence of dissipa¬ 
tive forces, that has been proved in Ref. [5] for the case of 
charged particle motion in an electromagnetic field, with 
the classical radiation reaction force taken into consider¬ 
ation. Particular examples of this effect were considered 
in the field of traveling Q and standing waves [7J. In 
the latter case, it is connected with stochastic heating 
(M3- However, due to the quantum nature, emission 
occurs randomly in the form of discrete photons that may 
randomize motion. As a consequence, the phase volume 
may be not compressed and even increases in spite of 


dissipation. Comparative analysis of an electron beam 
interaction with a counterpropagting laser pulse, taking 
into consideration radiative losses in classical and quan¬ 
tum form showed that at certain electron energies and 
laser field intensity, the phase space volume can be ex¬ 
panded during the interaction (llj . In this paper, we 
analyze the characteristics of electron motion in the field 
of a standing linearly polarized plane wave taking into 
account radiative losses in conformity with the quantum 
and classical descriptions. We consider the wave ampli¬ 
tudes at which an electron loses a small part of its energy 
in one act of photon emission but its motion can qualita¬ 
tively change due to discreteness of the photon emission. 
A comparative analysis of electron dynamics with and 
without radiation reaction force in classical and quan¬ 
tum forms is presented. 


ELECTRON MOTION IN THE CLASSICAL CASE 


For taking into account radiative losses in the classi¬ 
cal case, we add to the equation of motion the radiation 
reaction force in the Landau-Lifshitz form, keeping the 
most important term proportional to the squared Lorentz 
factor 12]: 
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where p, v, 7, — e, m are the momentum, velocity, 

Lorentz factor, charge and mass of the electron, corre¬ 
spondingly; c is the velocity of light; t is the time in the 
laboratory frame of reference; 5 = 2e 2 cj/(3rac 3 ); and uj 
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is the laser frequency. Without loss of generality, we as¬ 
sume that the electric field is directed along the Z axis 
and the magnetic field is directed along the X axis: 

E = E 0 cos (ut ± ipo) cos(ky)z, (3) 

B = E() sm(ut ± (fo) sin(%)x, (4) 

where the vacuum wave number k = u/c and cpo is the 
initial wave phase. The electron density is considered to 
be low enough, so that the effect of the particles influ¬ 
ence on each other can be neglected in comparison with 
the impact of the external field. Dimensionless variables 
Tj = ky , £ = kz, p = p/mc, r = ut, v = v/c, ao = 
eEo/(muc) are used for simplicity; we also assume that 
at the initial moment of time, the particle moves in the 
ZY plane. In this case, its motion always occurs in this 
plane: 

7-^r = -Pzdo sin(r + ip 0 ) sin (rj) - 5p y alQ{r, rj), (5) 

7= Pv a o sin ( T + <Po) sin (rj)- 
7 «o cos(r + ip 0 ) cos(r?) - 5p z a 2 0 Q(j , rj), ( 6 ) 

* = P>/1, ( 7 ) 

f = fc/7 - <8) 

where 0 = (y 2 — p 2 z ) cos 2 (r ± ipo) cos 2 ( 77 ) ± p 2 sin 2 (r ± 
^o)sin 2 ( 7 ?)-^-sin( 2 (r + ^o))sin(27?), p= ^l + z? 2 +p 2 , 

7 = yj\ ± p 2 . To analyze the electron motion with¬ 
out radiation reaction force it suffices to assume 0 = 0 . 
Note that, although further analysis is performed for the 
wavelength A = 0 .8pm, when a 1 , the system de¬ 
pends to a good accuracy only on one similarity param- 

eter ^3 = 4^ a 3 0 . 

In a general case, analytical solution of the equation 
of motion ([5]l8|) in the field of a plane linearly polarized 
standing wave hasn’t been found, even without radiation 
reaction, so the spectrum of characteristic Lyapunov ex¬ 
ponents is used to analyze electron motion with the help 
of a modified generalized Benettin algorithm 0,0. In 
order to find all possible modes of motion, it is neces¬ 
sary to consider all possible initial conditions. Because 
of the periodicity of the equations along the Y axis it is 
sufficient to set the initial position of the electrons along 
the axis within one wavelength. Taking into account the 
radiative reaction force, it is reasonable to assume that 
the absolute value of the transverse momentum (along 
E) during steady motion doesn’t exceed 2a, and the ab¬ 
solute value of the longitudinal momentum doesn’t ex¬ 
ceed a 2 / 2. However, significant computing resources are 
needed for good resolution of phase space by initial con¬ 
ditions. We confine ourselves to a single particle, which 
has zero momentum at the initial moment of time and is 
located between the electric and magnetic field nodes (in 


the calculations, the particle was set midway between the 
nodes). Thus, we haven’t obtained all possible modes of 
motion (for example, stable trajectory in the region of the 
electric field antinode at a relativistic initial momentum 
along the electric field 0), but such initail position and 
momentum allow the particle to sit on attractors, which 
are discussed in the work jjj. To reduce the influence 
of the initial conditions, the calculation of the Lyapunov 
characteristic exponents starts only after 30000 periods 
of laser field. As shown by calculations, this time is suf¬ 
ficient for the particle trajectory to sit on the attractor, 
if it exists, even with a ~ 1. Figure [T] shows the ob¬ 
tained spectra of Lyapunov characteristic exponents as a 
function of a. First, note that, as the dimension of sys- 



FIG. 1 . Lyapunov characteristic exponents of system ( 151171 ) 
with (a) and without (© = 0 ) (b) radiation reaction force. 
According to value/no CTapmuHCTBy Ai (red), A2 (black), A 3 
(blue), A4 (magenta). The sum of the exponents is X g (green). 


tern dSHH]) is 4, there exist chaotic attractors. Second, the 
results of numerical modeling are consistent with the the¬ 
orem of phase space contraction [H|], taking into account 
the classic radiation reaction: the sum of the Lyapunov 
exponents is negative, despite the fact that the maxi¬ 
mum Lyapunov exponent can be positive, the Lyapunov 
exponents equal zero. Consequently, system <1511511 has a 
chaotic or a regular attractor. Note that there is a range 
of amplitudes 35 < a < 138 where the attractor is reg¬ 
ular, which is conformed by the conclusion in [7]. The 
stochastic heating is suppressed in this amplitude range. 
Since the sum of Lyapunov exponents is negative, one 
of them is assumed to be zero, which improves accuracy 
of estimations of the other Lyapunov exponents. Average 
squared error of the exponents in this case did not exceed 
1 %. 

But without the radiation reaction, system ([5]t7J) may 
have attractors under certain initial conditions |8|, |16l |. 
However, the calculations havn’t given us definite infor¬ 
mation at what amplitudes attractors may arise. The 
exact value of the sum of the exponents X g could not 
be established because A^ 10 9 ± 0.01, i.e. the er¬ 
ror is much greater than the value itself. This occurs as 
the largest and the smallest absolute values of the Lya¬ 
punov exponents are opposite in sign | Ai ? 4 | ~ 0.2 ±0.002, 
and the second and third parameters are close to zero 
A 2j 3 ^ 10 -6 ± 10 -4 . However, the value of the sum of the 
Lyapunov exponents isn’t of great significance, most im¬ 
portant is that the senior indicator is positive and there 
is chaotic motion without radiation friction |7]. 
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Note that the spectrum of the characteristic Lyapunov 
exponents is calculated for a certain trajectory that orig¬ 
inates under certain initial conditions. However, accord¬ 
ing to the multiplicative ergodic theorem 0, this spec¬ 
trum can be attributed to the attractor as a whole. We 
have considered the attarctor which electron initially at 
rest reaches. Moreover if the initial momentum is po <C a, 
the particle goes to the same attractor [7]. We also note 
that access to the attractor, as follows from numerical 
calculations, does not depend on the initial phase of the 
wave. The initial phase affects the steady component 
of the pulse along the electric field (p z ), but due to ra¬ 
diation losses the steady component disappears in some 
time that will be estimated later. As a result, the average 
value of p z is equal to 0, and that’s why the initial phase 
has no influence on a possibility of reaching the attrac¬ 
tor. However, there are points in phase space starting at 
which an electron cann’t reach the attractor. The mag¬ 
netic field antinodes ys = 0.25A + 0.5nA (where n is an 
integer) are stable at small perturbations. In the non- 
relativistic case, the electron motion near this point can 
be described as damped oscillations in the ponderomo- 
tive potential due to radiation reaction force. It can be 
assessed in dimensionless units that \p\ <C 1 along with 
\y — Vb\ <C7 t/2 and a > 1 are sufficient for stability. The 
momentum can be assessed as p y ~ p z < \po\ + 2a\y — yB\, 
since a trajectory is rotated in a magnetic field with con¬ 
secutive increasing and decreasing of energy due to the 
electric field work. Therefore, a particle with zero ini¬ 
tial momentum is captured by the antinode of the mag¬ 
netic field, if 1 2/0 - Vb\ < 1/2 a. Due to the small size of 
this area around the antinode of the magnetic field, from 
which the particle will not go away, only a small frac¬ 
tion of the uniformly distributed particles is captured by 
this area. There is also an unstable point of equilibrium 
ys = 0.5nA, however, the arbitrarily small perturbations 
allow the particle to reach the attractor. Thus, only a 
small amount of particles does not reach the found at¬ 
tractor. In the case without radiation reaction, the ini¬ 
tial phase of the standing wave comes into play. For 
comparison of the results of the calculation of Lyapunov 
exponents with and without radiation reaction, the ini¬ 
tial phase is chosen to be 0, which corresponds to a zero 
steady component of the transversal momentum. With¬ 
out radiation reaction, the antinodes of the electric field 
are also unstable. The magnetic field antinodes are in 
indifferent equilibrium. Starting with an arbitrary po¬ 
sition outside a small region around the magnetic field 
antinode, the trajectories are qualitatively similar (for a 
long period of time, the trajectories fill the same space 
of phase projection p y p z y ). So, the Lyapunov exponent 
spectrum of a trajectory can be generalized to the others. 

Note that each characteristic parameter A^ corresponds 
to the basis vector i and along these vectors the trajectory 
perturbation increases or decreases as e Ait , depending on 
the sign of AAccording to numerical simulation, the 



FIG. 2 . Trajectory of an electron in the field of a linearly po¬ 
larized standing wave with amplitude a = 136 initially placed 
between electric and magnetic field nodes with zero momen¬ 
tum after 3000T and 7 and x as a function of time without 
radiation losses. 





FIG. 3. Trajectory of an electron like in Fig[2]but with clas¬ 
sical radiation reaction force. 

same basis vector (|9]) corresponds to the highest nonzero 
Lyapunov exponent with and without radiation reaction. 
The vector 77 is determined to an accuracy of 5% in the 
range of standing plane wave amplitudes 35 < a < 138. 
Disturbances along this vector decrease most slowly (the 
Lyapunov exponent is negative, the radiation reaction is 
taken into account) or increase most rapidly (the Lya¬ 
punov exponent is positive, the radiation reaction is ne¬ 
glected): 

r b = 0-ipy j ? Oit, Oiy) (9) 

Compare the dynamics of electrons in the field of a stand¬ 
ing plane wave with and without radiation reaction. We 
consider the wave amplitudes at which a regular attractor 
is possible assuming the radiation reaction to be taken 
into account. If the attractor is chaotic, no significant 
difference is observed in the motion of an electron with 
and without radiation reaction force. Moreover, the max¬ 
imum Lyapunov exponents are approximately the same 
in these two cases. For numerical analysis it is conve- 
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nient to consider the wave amplitude corresponding to 
the shortest time needed for the particle to approach the 
attractor; it occurs when a = 136. In this case, the 
largest non-zero Lyapunov exponent A 2 = —0.0067. Fig¬ 
ure 12131 shows the trajectories of the electron in the field 
of a plane wave with amplitude a = 136, with and with¬ 
out radiation reaction taken into account, respectively, 
after 3000 field periods. This time is sufficient for the 
particle motion to converge to the attractor. The parti¬ 
cle motion along the attractor is directed, the averaged 
transversal momentum is zero (there is no drift in the 
transversal direction), and the longitudinal momentum 
oscillates around ~ 150rac with an amplitude of ^ lOOrac 
and an averaged longitudinal velocity of 0.92c. Note that 
the averaged value of energy is < 7 amc 2 . Without 

radiation reaction, the electron motion is partially simi¬ 
lar to the electron motion with radiation reaction taken 
into account, but there is a possibility for the electron to 
be trapped by a magnetic field antinode. When trapped, 
it stays in this state for random time, drifting along the 
electric field, thus making the motion chaotic. This is 
the reason of directed motion loss. So, if we consider 
the electron bunch motion in the field of a standing wave 
with allowance for the radiation reaction, then we will 
have a M wave-like solution”’ (the electrons are halved 
and move in opposite directions along the Y axis). In 
this case, after a while there are no electrons in the ini¬ 
tial region. However, without radiation reaction, electron 
bunch spreading is observed (a M diffusion like solution”) 
due to the random change of electron motion direction. 
Let us analyze the electron bunch motion. 1000 elec¬ 
trons are initially uniformly distributed within the region 
—0.25A < y < 0.25A and — 0 . 1 A < z < 0 . 1 A. The initial 
momentum is zero. Figure 0]|6] shows the results of cal¬ 
culations. The quantum parameter x 1 d3,d3, i.e. 
the classical approach is applicable for description of the 
electron motion: 


X = -yj J (TOC 7 E + p x H) 2 - (p • E ) 2 (10) 

The energy spectrum of the electrons, with the radi¬ 
ation reaction taken into account, is discrete, and the 
spectrum line thickness is less than 0.2MeV. The num¬ 
ber of spectrum lines is determined by the time needed 
for the electron to approach the attractor, as well as by 
the fact that the electron motion is quasiperiodic and 
its period equals several periods of the standing wave 
(~ 14 in the case considered). The average energy of the 
electrons along the attractor is ^80MeV. The projection 
of the phase space portrait on the p y p z plane is a line 
that is the result of contraction of the phase space in 
FigHl In particular, this figure confirms that the major 
part of particles initially uniformly distributed within the 
—0.25A < y < 0.25A range converge to the same attrac¬ 
tor. In other words, for the majority of particles their 




FIG. 4. Elelctron energy distribution in the field of a plane 
wave with amplitude a = 136 (a), and phase space projection 
on the p y p z plane (b) with classical description of radiation. 



FIG. 5. Electron energy distribution in the field of a plane 
wave with amplitude a = 136 (a), and phase space projection 
on the PyPz plane (b) without radiation reaction. 


initial position is of no significance, which justifies our 
choice of the initial conditions for estimating Lyapunov 
exponents. The phase portrait without radiation reac¬ 
tion is ”blurred”, the spectrum becomes continuous and 
shifted to lower energies due to a random change of elec¬ 
tron motion direction, and the average energy decreases 
to 60 MeV (Fig. [5]). In the absence of radiation reaction 
force, several particles can reach higher energies. Part 
of the particles have 7 > 250, which is impossible if an 
electron moves along the attractor. 

The electron energy spectrum reflects the specific fea¬ 
tures of particle motion. There are divergent electron 
beams along the characteristic curves vt — y with radia¬ 
tion reaction. The electron distribution along the Z axis 
(~ 20 A) is determined only by the initial stage of ap¬ 
proaching to the attractor in Fig. [ 6 H), 0 and doesn’t 
change after that. Thus, the electron temperature tends 
to zero. Note that in ^ 1000T, the particles escape from 
the area where they were initially located. Figure EH) 
confirms that the position of the leading edge of the elec¬ 
trons corresponds to 0.92c£. If the radiative losses are 
neglected, the diffusive spreading of the electron distri¬ 
bution along the Y and Z axes is observed due to particle 
trapping by the antinode of the magnetic field and chang¬ 
ing direction of motion (Figj 6 H> ),[!]). In this case, the par¬ 
ticle motion can be described by the diffusion equation 
for the linear density of the particles: 
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FIG. 6. Time dependence of electron distribution along Y 
and Z axes for a = 136 with (a) and without (b) radiation 
reaction force. 
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FIG. 7. Electron distribution along Y and Z axes after 2000 
(red dashed line) and 5000 (black dashed line) wave periods. 



FIG. 8. Diffusion of electrons along Y (a) and Z (b) 
axes. Black and red lines correspond to time t=2000T and 
t=5000T. Chain line corresponds to numerical calculation of 
1000 trajectories of electrons initially distributed in the A/2 
range along Y and Z axes. Solid line corresponds to solution 

m 


diffusion coefficient D y ^ 120, and the transverse one 
D±~ 1. 

Comparison of the electron motion with and without 
radiation reaction points that the qualitative differences 
in the motion occur only after the particles approach the 
attractor, and it takes hundreds of wave periods and con¬ 
tinuous action of dissipative force in time. However, an 
electron emits a portion of its energy, quanta, and each 
act of emission "‘knocks’”(or pushes) an electron, thus 
leading to lengthening of the time of convergence to the 
attractor or disappearance of the attractor at all. To un¬ 
derstand the influence of the discreteness of photon emis¬ 
sion on electron motion, we can estimate how often the 
electron emits and what portion of its energy turns into 
radiation. For this purpose it is necessary to estimate 
the parameter %. If the electromagnetic field is specified 
by (|3l4lh % can be written as 


where r is implied to be y or z. As can be seen from 
Fig. [ 6 l the diffusion along the Y axis is much faster, so it 
is possible to consider the diffusion in each direction in¬ 
dependently. Since the electrons are initially within the 
half wavelength range and their following motion is con¬ 
sidered during 5000 periods of the standing wave, during 
this time a typical scale of electron distribution becomes 
A. Then, it is possible to use the fundamental solution, 
assuming the diffusion coefficient to be constant and the 
initial distribution to be a delta function multiplied by 
the number of particles (fT 2 |) : 


N er 


No 

— . : exp 

v/47rD|| j± t 



( 12 ) 


where N eyjZ is the linear electron density along the Y,Z 
axes. The evolution of the electron density profile along 
the Y and Z axes is shown in Fig. [ 8 j Calculations show 
that, if the wave amplitude a = 136, then the longitudinal 


X = — 2 a \/(pl + !) cos 2 (r + <p 0 ) cos 2 (??) 
mc z v y 

+p 2 sin 2 (r + (fo) sin 2 ( 7 ) 

- 0 . 5^7 sin( 2 ^) sin( 2 (r + <p 0 )) (13) 

Rough estimation is % ^ ^3 a 2 ~ 0-03, which demon¬ 
strates a good fit with the value of x along the tra¬ 
jectory with and without classical radiation reaction 
(Fig l2l3D . In the ultrarelativistic case, an average pho¬ 
ton energy can be estimated as e p h/mc 2 ~ 0 . 457 % ~ 
a 3 huj/mc 2 rsj 3.3 0 , and the classic radiation power 
Prad = holds valid in the case % << 1 . 

Hence, the average number of photons during the period 
can be assessed to be N ph « ^^ 77 “ ~ ~ 

and part of electron energy lost in one act of photon 
emission is %/ 2 . Note that, as follows from calculations, 
each of 1000 electrons in the wave field with amplitude 
a = 136 during 5000 field periods emits N p h = 6 on the 
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average. Rough estimations are in a good agreement with 
numerical results. Let us determine an average number 
of photons emitted by an electron over the wave period 
according to the probability of photon emission per unit 
time in the quantum case [ 2 1 ] : 


role. Note that the system of equations ([5]|8]) when a 1 
is determined by one similarity parameter 5a 3 oc c da 3 . 
Thus, in contrast to the classical case, there is no reg¬ 
ular attractor and suppression of stochastic heating, if 
(0.0006/T) 1 / 3 < a < (0.037/T) 1 / 3 and cj < 0.05c/r e or 
A > 100 r e « 0 . 001 A 


W = 


1 e 2 me 2 
3\/3?r he fry 


/ 


5 u 2 + 7u + 5 

(1 + m ) 3 


K 2 /^u/2>x)du. 


(14) 

When x <C 1, the expression d reduces to W = 
2 ^ 377 ^hj X • This an estimate of a photon num¬ 
ber N p h = ~ 9, and each act of photon emission 

takes away about 0.02 part of electron’s energy. Thus, 
(27r£a 3 ) -1 ^ 5 wave periods are sufficient for an elec¬ 
tron to lose a steady component of the momentum along 
the electric field. It may seem that the classical descrip¬ 
tion can correctly describe electron motion due to a small 
fraction of electron energy converted to photon energy, 
but discreteness of emission plays an important role. On 
the one hand, the energy loss must lead to formation 
of an attractor. On the other hand, randomness leads 
to divergence of the particles in phase space. So, exis¬ 
tence of an attractor is determined by the competition 
of these two factors. Estimations show that the parti¬ 
cles move about 0.17 wave periods between the acts of 
photon emission without radiation reaction, and in this 
case the maximal Lyapunov exponent is Ai = 0 . 2 . This 
means that the particles which are close in phase space 
diverge in ~ T/(27rAi) ~ 0 . 8 T. The time of particle con¬ 
vergence due to the radiation reaction may be estimated 
to be ~ T/(2 ttA 2 ) ~ 24T. The difference of these times 
destroys the attractor, so it is necessary to consider ra¬ 
diation as a quantum mechanism. 

The importance of a quantum mechanism may also 
have a different explanation. Suppose that there is a reg¬ 
ular attractor with allowance for quantum description of 
photon emission. Then, the attractor will be like in the 
classical case neglecting discreteness of radiation, since if 
X <C 1 , then radiation losses are the same in the classical 
and quantum cases. Hence, the value of the characteris¬ 
tic Lyapunov exponents must be the same too. In this 
case, as was shown earlier, the largest non-zero Lyapunov 
exponent is |A m ##| <C 0 . 1 , which corresponds to the ba¬ 
sis unit vector (|9]) . Perturbation along this vector during 
the time period between the radiation events A T ra d de¬ 
creases as e _ l AmRi? l ATrad . However, taking into account 
the discrete nature of radiation, it follows that the par¬ 
ticle during A r ra d is affected only by the Lorentz force, 
and in this case the largest non-zero Lyapunov exponent 
A md ~ 0.2 \XmRR\i and the perturbations grow along 
the vector (J9]) as e XrndArrad . The result is that the pertur¬ 
bations grow as e ( Amd+AmRR ) Ar , since A m d + A m RR > 0, 
and a regular attractor cannot exist. Consequently, the 
discreteness of radiation plays a qualitatively important 


ELECTRON MOTION TAKING INTO ACCOUNT 
DISCRETENESS OF RADIATION 

Let us use the tried-and-true approach for the quan¬ 
tum description of radiation and classical description of 
electron motion [H]. Electron motion is described by 
equations 11121 but without radiation reaction. Each par¬ 
ticle should propagate a randomly assigned optical length 
in the electromagnetic field before photon emission. If 
the random number rq G (0,1), then the optical path 
L op t = log—Once fWdt = L opt , then the pro¬ 
cedure of photon emission is started. The choice of 
the photon energy is carried out with the help of the 
Monte Carlo Rejection sampling method j23|. To do 
so, the auxiliary function M(rj) is selected, such that 
M(rj) • C(x) > dW(Xi v) /dv where 77 = e p h/{mc 2 7 ) is 
the fraction of the electron energy converted to the pho¬ 
ton energy: 


M( V ) = 


' 0.325 
^2/3 5 

0.18 


77 < 0.87 

(i-„)i/a > 0.87 < rj <1 


(15) 


C(x) = 1.6 


e 2 me 2 2/3 

he hwx 


(16) 


Here, x is calculated according to the momentum and 
strength of magnetic and electric fields in the constant- 
field approximation; dW{x^ r f)/dr] is defined as follows 

El: 


dr] 


1 e 2 me 2 
7 t >/3 he hwx 


/ 

2t7 

_3(1-77)x 


^ 5/3 (y)dy+ 


1 — ?7 


^2/3 ( 


2 V 


3(1 - V)X 


(17) 


The second random number r 2 G (0,1) allows deter¬ 
mining rj(x, r 2 ) from the equation M(rj) = r 2 . Us¬ 
ing the third random number r% G (0,1) it is possi¬ 
ble to determine whether the photon is emitted rs < 
dW (x, 17 )/ dr] / (C (x) * M(rj))\ otherwise, it is necessary 
to choose a random number r 2 again. A properly cho¬ 
sen function M(rj) is as close as possible to the func¬ 
tion dW(x , v)/dr]/C(x) : which reduces the number of at¬ 
tempts to emit a photon. Also, it is easy to determine 
77(x, ^ 2 ), he. to find the inverse function. 
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Let’s consider a typical trajectory of an electron 
using the quantum description of radiation (Figj9j). 
Trajectories in cases of quantum description and without 
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FIG. 9. Electron trajectory as in Fig[2] but with quantum 
description of photon emission. 


(FigH2J). Calculations demonstrate that the diffusion co¬ 
efficients are smaller than in the case without radiation 
friction, as the electron changes its direction more of¬ 
ten. Thus, the diffusion decreases and the electron bunch 
"spreads” slower. 



FIG. 11. Electron distribution along Y and Z axes as a func¬ 
tion of time for a — 136 using quantum description of photon 
emission. 


radiation losses are similar, as is evident from compari¬ 
son of Figs. [9j [2j but there is one qualitative difference. 
Without radiation reaction, there is an integral of mo¬ 
tion p z — A(y,t) = const, because the equation of mo¬ 
tion does not depend on the coordinate z. Therefore, the 
motion along the z coordinate is different for different 
initial phases of a standing wave, the electrons can have 
a steady component of momentum along the Z axis, but 
this component vanishes both with allowance for classical 
and quantum radiation losses. 



FIG. 10. Electron energy distribution in the field of a plane 
wave with amplitude a = 136 (a), and phase space projection 
on the p y pz plane (b) with quantum description of photon 
emission. 


The longitudinal diffusion coefficient is D\\ ~ 100, and 
the transversal diffusion coefficient D± ^0.6. In order 
to determine why the diffusion is like this, note that the 
motion of an electron in a plane wave with quantum de¬ 
scription of radiative losses is chaotic. It is found that 
the electron retains the character of the motion each half 
period, and the character of motion in the next half pe¬ 
riod is determined by the present one. Such motion can 
be regarded to be a Markov chain. 




Let us now examine electron bunch motion under the 
initial conditions same as in the classical case. Qualita¬ 
tively, the motion of the electron bunch coincides with 
the case without radiation friction, which confirms the 
conclusion of the importance of discreteness and that 
the motion is generally determined by the Lorentz force 
rather than the radiation reaction (Fig jlOIlT]) . Figure 
nn shows that part of particles may have energy greater 
than 250mc 2 . The spectrum has an analogous feature, 
when the radiation reaction is not taken into consider¬ 
ation. This is due to the discrete nature of radiation, 
because particle does not undergo radiation damping be¬ 
tween the acts of emission. The electron bunch motion 
can also be described using the fundamental solution [T2] 


FIG. 12. Diffusion of electrons along Y (a) and Z (b) axes 
at quantum description of photon emission. Black and red 
lines correspond to time t=2000T and t=5000T. Chain line 
corresponds to numerical calculation of 1000 trajectories of 
electrons distributed initially in A range along Y and Z axes. 
Solid line corresponds to solution GUi 


MARKOV CHAIN 

The process described by the Markov chain is defined 
by a transition matrix which sets the probabilities of the 
states of the n th step from the known state of the n-l th 
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step. We use the formalism of Markov chains to deter¬ 
mine the width of the electron distribution along the Y 
and Z axes. There are 3 states of motion along the Y axis 
(all the states are shown in Fig. fl3l) . 1) A particle shifts 
by 0.45A in the direction opposite to the Y axis, while 
shifting by T0.22A along the Z axis. The probability of 
changing the direction with respect to Z without chang¬ 
ing the motion along Y is 0.35. The probability of the 
transition to the second state is 0.023. 2) On the average, 
the particle propagating T0.25A along Z does not move 
along Y. The probability of changing the direction of mo¬ 
tion along the Z axis is 0.4. The probability of transition 
to state 3) or 1) is 0.15. State 3) is the same as state 
1), but the electron moves along the Y axis. Numeri¬ 
cal simulations allow determining the probability of the 
transition between the states (indicated above the arrows 
in Fig. USD, as well as determining the states. Assume 

/s?\ 

that S n = s? is the state vector of electron motion 
\SS ) 

along the Y axis corresponding to the n th step. S% and 
S 3 denote the probability of the motion along the Y axis 
and in the opposite direction with respect to the Y axis 
correspondingly, S 2 denotes the probability of rest along 
the Y axis. Then the transition matrix Pij from state i 
to j is written in the form 

/0.977 0.023 0 \ 

Pij = 0.15 0.7 0.15 (18) 

\ 0 0.023 0.977/ 

which means that P n =(P T ) n S°. Let’s calculate 



FIG. 13. Schematic electron trajectory as Markov chain 


the mathematical expectation /x and the variance 
D of electron deviation along the Y axis in n steps 

n n 

AY E = = E(—0.45A,0A,0.45A)5L The mean 

i=1 i=1 

value of the sum is the sum of the mean values, i.e. 

n 

/jlSYz = E 0.45A, 0A, 0.45A)5T Initially, the parti¬ 
al 

cles are distributed uniformly in states 1), 2), and 3) 
according to numerical simulations of electron bunch 
n f 1/3 \ 

dynamics, i.e. S {) = 1/3 . /aSY^ = 0 because of the 

V 1/3/ 

equivalent motion in both directions along the Y axis 
and because of the symmetry of matrix (P T ) n and 
the vector of initial states. Thus, as evidenced by Fig. 
na on the average, motion occurs without drift of the 
center of mass. Let’s calculate the variance PAYe = 


J EE riW) = ./E M*? + 2 E E KYiYj)- Note 

y*=ii=i y »=i i=ij>i 

that the events become weakly correlated after a suf¬ 
ficiently large number of steps and the state vector at 
the k -th step ( P T ) k S° = S av . In order to find S av , it 
is necessary to solve the equation P T Pj v = Pj v in con¬ 
junction with the condition that the sum of the values in 
the rows of equilibrium transition matrix P av (a matrix 
with three identical columns) is equal to 1. Because of 
the symmetry of motion in both directions along the 
Y axis, the equilibrium transition matrix has the form 
f a b a^ 

2a + b = 1, and 0.977a + 0.156 = a. 


P — 

1 av — 


aba 
< a b a j 


It 


Hence, a = 0.4644, b = 0.0712, and S av = 
took N eq « 500 steps for the equilibrium state to be 
formed. The electron bunch dynamics is observed during 
5000 wave periods, i.e. N st = 10000 steps. The first 
term in the expression of variance can be calculated as 

n 

E ME 2 = = 0.188nA 2 , where Yi is the displace- 


i =1 

ment during half the wave period along the Y axis. Let’s 
calculate the second term. Assume that S l is the state 
vector at the i th step. For finding the correlation of 
states i and j, it is necessary to know the state vector S J 
using the conditional probability, if we know the value of 
random variable Y {. These probabilities are expressed in 
the form of elements of the matrix (P/^)- 7- *, where the 
I th column indicates which states will be at the ] th step. 

n n 

If the state k is at the \ th step, then E E MEY)) = 

i=l j>i 



-°of A ) T (p-r f 

0.45A / V ' V 


-0.45A 0 0 

0 0 0 
0 0 0.45A 


) (. p T ys°. 


The approximation ( P T ) k S° = S av is used because 
N st > N eq , hence A = ( °’o 45A 0 0 ) (P T YS° = 

( b ). Note that (P T YA tends to zero for large 
V 0 . 21 A ) 

i, since components 1 and 3 of the vector A are 
similar in magnitude and opposite in sign, in other 
words the dependence on initial conditions disappears 
and the motions along the Y axis are equal in both 
directions. Thus, for good accuracy it is sufficient 
to take into account a certain quantity of terms 
(. P T ) l A . It suffices to take into account only 500 
terms in view of the fact that the equilibrium state 
is achieved in 500 steps. Finally, to a good accuracy 


n n 500 / _ n 4 n x T 

E E K Y iYj) = «E( oa ) (P T YA = 7 . 9 n\ 2 . 

i=l j>i j=1 ' 0.45A / 

Generally, the variance is due to correlations and 
PA Ye = Ay/nX. Based on the root function of the 
number of steps (twice the number of periods) it becomes 
obvious that M spreading" of an electron bunch occurs at 
large n with constant dimensionless diffusion coefficient 
P|| = 4 2 • 2tt ~ 100. The value of the variance is in a good 
agreement with the characteristic width of the electron 
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posed for describing electron kinetics in this case using 
Markov chains. The electron distribution obtained by 
this method is in a good agreement with the results of 
direct calculations. 

The research is partly supported by the grant (the 
agreement of August 27, 2013 N02.B.49.21.0003 between 
The Ministry of Education and Science of the Russian 
Federation and Lobachevsky State University of Nizhni- 
Novgorod) and RFBR (grant NT4-02-31495 mol_a). 


FIG. 14. Electron distribution over 2000 wave periods (red 
line) and over 5000 wave periods (black line) along Y a) and 
Z b) axes as a result of Markov chain process (solid line) and 
as a result of numerical solution of motion equation (dashed 
line) taking into account quantum description of radiation. 

distribution along the Y axis in Fig. [12] and with the 
diffusion coefficient. Figure [IT] shows the distribution 
of electrons along the Y and Z axes derived from the 
Markov process according to Fig. [13] in comparison with 
the numerical simulation of electron motion taking into 
account the quantum description of radiation. 

Similar analysis can be made with the electron motion 
along the Z axis. In this case, the state vector consists 
of four elements. The result is the following: DAZ ^ = 
0.3^/nX and D± = 0.3 2 • 27t ~ 0.6, which agrees well with 
the results presented in Fig. [12] 

It should be noted that the consideration of the elec¬ 
tron motion as a process described by a Markov chain 
allows us to understand the transition to the normal ra¬ 
diative trapping 0 . 0 - The state of Yi = 0A becomes 
absorbing, i.e. in terms of the transition matrix we can 
say that P 22 —> 1, P 21 and P 23 tend to zero with increas¬ 
ing a. 


CONCLUSION 

The electron motion in a plane standing wave is con¬ 
sidered within the framework the quantum and classical 
description of the radiation reaction. It is shown that 
in the case of small y, the classical description of pho¬ 
ton emission can lead to qualitatively wrong electron dy¬ 
namics despite the small radiative losses. The Lyapunov 
characteristic exponents obtained in this case show that 
there is a regular attractor in the phase space, whereas 
without radiative losses the motion is chaotic. In the 
quantum case, the character of motion changes drasti¬ 
cally due to the discrete nature of photon emission, which 
randomizes particle motion. This leads to the additional 
divergence in the phase space, which completely destroys 
formation of a regular attractor. A new method is pro- 
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